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Abstract 



We consider the phenomenology of GMSB models where the lighter stau is 
the next to lightest supersymmetric particle. In this situation, the dominant 
signals for supersymmetry are events with two or three high T leptons 
accompanied by large missing transverse energy. We find that the inclusive 
two r-jets signature could be observable at the Tevatron's Run II, while the 
inclusive three r-jets signature could be important at Run III. 
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I. INTRODUCTION 



The phenomenology of gauge mediated supersymmetry breaking (GMSB) models have 
been the subject of much interest lately [0-01 . They provide an alternative to the usually 
considered case of the soft terms involving the low energy fields induced by gravity. In 
GMSB models the supersymmetry (SUSY) breaking is communicated to the visible sector 
by gauge fields. The scale at which this occurs is usually taken to be around 10 5 GeV. 

The sparticle spectrum in GMSB models has some significant differences from the usual 
gravity mediated SUSY breaking models. In GMSB models, the gravitino is the lightest 
supersymmetric particle (LSP). The next to lightest supersymmetric particle (NLSP) is 
typically the lightest neutralino or the lighter stau. Most phenomenological studies and 
experimental searches that have been performed in the context of GMSB have taken the 
lightest neutralino to be the NLSP. The lightest neutralino then decays by Xi ~^ jG. If this 
decay takes place within the detector, the signal involves high p? photons accompanied by 
large Ifix [[§• For much of the parameter space, however, the lighter of the two scalar staus 
is the NLSP. In this case, the decays of SUSY particles produce the f\ which subsequently 
decays to a r and the gravitino. If the f\ decays occur within the detector, signatures for 
SUSY production will then generally include r leptons from the f\ decays and JpT due to 
the stable gravitinos and neutrinos leaving the detector. 

It was proposed that GMSB models where the t\ is the NLSP can lead to unusual 
and distinguishing signatures for SUSY production at colliders []5J. At the Tevatron, these 
signals arise from chargino pair production (xtXi) an d from the production of the chargino 
with the next to lightest neutralino (XiX:])- Slepton pair production (ifli where I = e, 
/i, or r) can also be significant. The subsequent decays typically involve high p T r leptons 
and substantial $t- in this paper we analyze the signals for such SUSY production at the 
Tevatron. In particular we consider the parameters for the n = 2 stau NLSP model line of 
the Gauge Mediation/Low Scale SUSY Breaking working group of the Physics at Run II 
Supersymmetry/Higgs workshop. 

II. MASS SPECTRUM AND PRODUCTION MECHANISMS 

Since the observed signal depends on the masses of the sparticles, we first begin by 
describing the model and the corresponding mass spectrum. In our model, the messenger 
sector consists of some number of multiplets that are 5 + 5 representations of SU(5). They 
couple to a chiral superfield S in the hidden sector whose scalar component has a vacuum 
expectation value (VEV) (s) and whose auxiliary component has a VEV (F s ). By imposing 
the requirement that the electroweak (EW) symmetry is broken radiatively, the particle 
spectrum and the mixing angles depend on five parameters: M, A, n, tan j3 and the sign of 
fi. M is the messenger scale. A is equal to (F s ) / (s) and is related to the SUSY breaking scale. 
The parameter n is dictated by the choice of the vector-like messenger sector and can take 
the values 1, 2, 3, or 4 to satisfy the perturbative unification constraint. The definition of 
tan (3 is taken as tan (3 = u 2 /i>i where V2 is the VEV for the up-type (H u ) Higgs doublet and 
v\ is the VEV for the down- type (Hd) Higgs doublet. The parameter /x is the coefficient in 
the bilinear term, [iH u Hd, in the superpotential. Constraints coming from b — > 57 strongly 
favor negative values for \x in our convention || and, in the cases considered in this work, 
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fj, is taken to be negative. Demanding that the EW symmetry be broken radiatively fixes 
the magnitude of \x and the parameter B (from the Bu,H u Hd term in the scalar potential) 
in terms of the other parameters of the theory. 

The soft SUSY breaking gaugino and scalar masses at the messenger scale are given by 

m 

M i (M)=n,(A)^M A (1) 

and 

™ 2 (M)=2n/(A)f>a (^)'a 2 (2) 

where the an are the three SM gauge couplings and ki — 1, 1 and 3/5 for SU(3), SU(2), 
and U(l), respectively. The Cj are zero for gauge singlets and are 4/3, 3/4 and (Y/2) 2 for 
the fundamental representations of SU(3), SU(2) and U(l), respectively (with Y given by 
Q = Is + Y/2). g(x) and f(x) are messenger scale threshold functions. We calculate the 
sparticle masses at the scale M using Eqs. (|1|) and (0) and run these to the electroweak scale 
using the appropriate renormalization group equations ||. 

The decay chains and hence the signatures for the events depend on the particles initially 
produced as well as the hierarchy of the masses. Since SUSY breaking is communicated to the 
visible sector by gauge interactions, the mass differences between the superparticles depend 
on the their gauge interactions. This creates a hierarchy in mass between electroweak and 
strongly interacting sparticles. Eq. (P shows that the gluino is more massive than the EW 
charginos and neutralinos, while Eq. (0) shows that squarks are considerably more massive 
than sleptons. Given this hierarchy of sparticle masses and the current lower bounds on 
squark and gluino masses, the production of strongly interacting sparticles is probably not 
a viable search modes for SUSY at the Tevatron Run II. A more likely mechanism for 
producing SUSY particles is via EW gaugino production. At the Tevatron, chargino pair 
{XiXi) production takes place through s-channel Z and 7 exchange and X2X1 production 
is through s-channel W exchange. Squark exchange via the t-channel also contributes to 
both processes, but the contributions are expected to be negligible since the squark masses 
are large in GMSB models. The production of x\ xt 1S suppressed due to the smallness of 
the coupling involved. In regions of the parameter space where the production of charginos 
and neutralinos is kinematically suppressed, the pair production of sleptons (fifi, fixfii and 
ei&i) can be important. Their production occurs through s-channel Z and 7 exchange. 

Given the hierarchy of sparticle masses in GMSB models, there are roughly four possible 
cases to consider for SUSY production at the Tevatron: 

Case 1: nip > M x o > M x ± > > M x o > m fl 

Case 2: M x o > M x ± > rxiy > M x o > > m fl 

Case 3: M x o > M x ± > rriy > m e - 1)(il > M x o > m fl 

Case 4: niy > M x o > M x ± > M x o > m e - li(il > m fl 
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The three sneutrino masses are nearly the same. The lighter of the selectrons and smuons 
are essentially right handed and have the same mass. Also, for all the parameter points we 
considered, Xx an d x\ are nearly degenerate in mass. 

The possible final state configurations at the Tevatron depend on the sparticle spectrum, 
but they will have certain aspects in common. Since the f\ is the NLSP, the various possible 
decay modes will (usually) produce at least two r leptons arising from the decays of the fx a. 
In addition, there can also be large Ifir due to the stable gravitinos and neutrinos escaping 
detection. 

A special situation in cases 2 and 4 arises when the lighter selectron and the lighter 
smuon are nearly degenerate in mass to the lighter stau. When the mass difference between 
the slectron and the stau is less than twice the mass of the r, essentially the only decay 
mode for the selectron is e — > eG. The lighter smuon likewise decays via jl — > \iG. This 
situation is refered to as the "co-NLSP" case ||. 



III. ANALYSIS AND RESULTS 

We now give a detailed analysis of the possible Tevatron signatures for SUSY production 
in the context of GMSB models where the lighter stau is the NLSP and decays promptly 
within the detector. This analysis is performed in the context of the Main Injector (MI) 
and TeV33 upgrades of the Tevatron collider. The center of mass energy is taken to be 
yfs = 2TeV and the integrated luminosity is taken to be 2fb _1 for the MI upgrade and 
30 ftr 1 for the TeV33 upgrade 0. 

In performing this analysis, the cuts employed are that final state charged leptons must 
have pr > 10 GeV and a pseudorapidity, r] = — ln(tan |) where 9 is the polar angle with 
respect to the proton beam direction, of magnitude less than 1. Jets must have E? > 10 GeV 



and \rj\ < 2. In addition, hadronic final states within a cone size of AR = y (A0) 2 + (Ar/) 2 = 
0.4 are merged to a single jet. Leptons within this cone radius of a jet are discounted. For 
a r-jet to be counted as such, it must have \r]\ < 1. The most energetic r jet is required to 
have E T > 20 GeV. In addition, a missing transverse energy cut of Ifi? > 30 GeV is imposed. 

We consider each mass case in turn. In our analysis, we restrict ourselves to those regions 
of the parameter space where the fx decays promptly to a r and a gravitino. The parameter 
space is also restricted to those regions where ^ 70 GeV. Results from LEP-2 set a 



lower bound of m fl > 72 GeV W2 



IV. STAU NLSP LINE WITH N = 2 

In this section we do the analysis for points along the line defined by the parameter 
values M/A = 3, n = 2 and tan/3 = 15. We vary A from 35TeV to 85TeV. The masses 
for the sparticles that are of interest here are given in Fig. |l|. Note that the sneutrino mass 
is always above that of the lightest chargino and the lightest two neutralinos. Thus the 
sneutrinos do not figure into the decay chains of the major SUSY production mechanisms 
here. Note that the lightest neutralino is below the selectron/smuon mass at the lower end 
of the A scale (A < 45TeV). Thus for A < 45TeV the mass spectrum is of type 1. For 
A > 45 TeV, the mass spectrum is of type 4. In this region of A, there are more decay modes 
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for the various particles due to the increasing masses of all the sparticles as well as the shift 
in the position of the lightest neutralino in the mass hierarchy. 

The cross sections for these parameters are given in Fig. ^|. From the figure, the cross 
sections for XiXi an d X^xt production dominate for the region where A is below 65 TeV. 
As A increases, the masses of the gauginos increase significantly and hence the cross section 
falls off. For A > 70 TeV, the production rates for fi ff~ and e\e{ are dominant, but the 
cross sections tend to be rather low. 

The signature depends on the allowed decay modes of the sparticles and their branching 
ratios. The branching ratios for the sparticles of interest are given in Table |. Since f\ is the 
NLSP, it decays via f\ — > r G. The decays of the selectron and smuon depend strongly on 
the value of A. For A below ~ 45 TeV, the lightest neutralino has a mass below that of fx\ 
and t\. As a consequence of this, the main decay mode of lightest smuon is Ji\ — > Xi A* an d 
the main decay mode of the lightest selectron is e\ — > Xi e - F° r higher values of A, however, 
the lightest neutralino mass increases above that of e% and fi\. Then the only available 
two-body decay mode for the lightest smuon is p,\ — > (iG with the corresponding decay 
mode for the lightest selectron. Given the smallness of the coupling involved, though, the 
possibility exists that some three-body decays could be important. Indeed, the neutralino 
mediated decays /ij" — ► e~r~f + and jl^ — ► e~r + f~ are the important decay modes |TT[ for 
these higher values of A. 

Since the lightest neutralino tends to be one of the least massive sparticles, its only 
decay modes are Xi an d the decays to [i\ and e\ if Xi is greater in mass than those 

sparticles. Since the lightest neutralino is less massive than the selectrons and smuons for 
A, the only decay mode is Xi ~ * ^i r - As A increases above 45 TeV, jl\ and e\ become 
increasingly important, although Xi ~^ remains the dominant decay mode. 

Since the lightest chargino is mostly wino, it couples mainly to "left-handed" sfermions. 
Thus the decay mode Xi ~ > ^i u t is typically important due to the significant mixing of the 
left-handed and right-handed staus and the lower mass of the f\. This decay mode is, in fact, 
essentially the only decay mode for low values of A with the parameters considered here. 
Thus with the subsequent decay f\ — > tG, there are typically two r leptons produced in 
XiXi production at these values of A. As A increases, however, the decay mode xt ~^ XiW 
becomes available and becomes the dominant decay mode as A increases above 60 TeV. With 
the two r leptons that can be expected from the lightest neutralino decay and the W — > tu t 
decay, we can expect up to six r leptons from XiXi production at these larger values of A. 

There are many decay modes for the second lightest neutralino as Table | shows. For 
low values of A, the dominant decay mode is xt ~> a ^ 50 - 60%. The decays to the other 
sleptons are also important at 15 - 20%. Thus xfxt production produces three r leptons: 
two from the slepton decays of the neutralino and one from the decay Xi ~^ T\ v t followed 
by t\ — » t G. As A increases above 55 TeV, The decay x\ ~^ Xi n ) where h is the Higgs 
boson, rapidly becomes the dominant decay mode. The decay X2 ~ * Xi % is a lso present, 
but of relatively little importance. 

Given the cuts that we place on the r-jets, the question arises as to how high we can 
expect the Et of the r-jets to be. Fig. Ogives the Et distribution of the highest Et r-jet for 
A = 35 TeV. The pseudorapidity cut of \r]\ < 1 on r-jets has been imposed in Fig. 0(b). The 
peak in the distribution occurs at about 25 GeV with a broad tail that reaches out beyond 
120 GeV. Thus the leading r-jets are relatively hard and many will pass the transverse 
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energy cut of Et > 20 GeV. The next to highest Et r-jet is significantly different as Fig. |] 
shows. Here the distribution peaks at a lower value of about 15 GeV and hardly extends at 
all beyond 80 GeV. Due to the softness of the secondary r-jets, many of the r-jets will tend 
to be eliminated by the cuts. 

Also of interest is the fa distribution. With energetic and stable gravitinos and neutrinos 
produced in the decays, it is expected that large missing transverse energy could be an 
important part of the signal. Since the missing transverse energy is calculated from what is 
observed, however, the question arises as to whether significant cancellation occurs due to 
the many decay products. Fig. |5| gives the fa distribution for the case where A = 35 TeV. 
The figure demonstrates that the fa distribution is indeed broad with a tail reaching out 
beyond 120 GeV. The peak before cuts occurs at about 35 GeV and the peak still occurs 
at about 35 GeV when Et/pt and pseudorapidity cuts are applied to the various particles. 
Thus a 30 GeV cut should not be too restrictive. As A is increased, the fa distribution gets 
harder since the gaugino masses get larger as A is increased. 

We now consider the specifics of the various final state possibilities. Table |J give the 
inclusive branching ratios for different number of r-jets for A = 35 TeV. As indicated above, 
this example always produces two r leptons in chargino pair production. Before cuts the 
inclusive branching ratio for the 2 r-jet mode in chargino pair production is 42%, while 
the 1 r-jet mode in chargino pair production is 45.6%. After the cuts specified above, the 
branching ratios are cut down rather substantially. The one r-jet BR becomes 25.8% and 
the two r-jet BR is 10.8%. The situation changes as A increases. This is demonstrated 
in Tables [III] and [V| which are for A = 50 and 70 TeV, respectively. We see that there 
is now the possibility for many more r-jets. This is due to the appearance of the decay 
mode xt ~ > Xi W. With the decay Xi ~^ followed by fi — > tG along with the decay 
W — > tu t , there is now the possibility for producing up to six r-jets. At A = 50 TeV, the 
BR for three r-jets is 16.3%. This is cut down substantially after cuts, however, and the BR 
becomes only 1%. A = 70 TeV shows similar results, although at this point this point the 
cross section for XiXi production is low. 

Considering X2 Xi production, we recall that at A < 45 TeV, the decay modes of X2 are 
X2 T\T, X2 ~ > A'i/- 4 an d X2 ~~ * ^i e - With the subsequent decays f\ — > r G and fi,i — > [Lrf\ 
along with the corresponding decay of the selectron, the number of r leptons from X2 decay 
is two. With the one r lepton from the chargino, we have up to three r leptons from X2X1 
production at these values of A. We see from Table [II] that for A = 35 TeV, the branching 
ratios for inclusive production of 3 r-jets is 27.2% before cuts, while for 1 and 2 r-jets 
the branching ratios are 24% and 44.3%, repsectively. These branching ratios are cut back 
considerably once the cuts are included. The 3 r-jet BR in particular is reduced to only 
2.6%. As A increases, the potential exists to create more than three r-jets due to the new 
decay modes for xt an d The branching ratios for more than three r-jets tend to be 
small after cuts, however, as Tables [TTI| and [TV] demonstrate. 

Slepton production tends to be rather simple since there are relatively few decay modes 
available to the sleptons. This is especially true for the production of the lighter stau as its 
only significant decay mode is f\ — > r G. Thus for f^rf production up to two r-jets are 
possible. Before cuts, the 2 r-jet BR is 42% and the 1 r-jet BR is 45.6% for f^rf production 
independent of the values of A considered. After cuts, these drop down to 10 - 17% and 30 
- 40%, respectively, as shown in Tables [TJ, |TJ and [TV]. The three-body decay modes of the 



6 



lightest selectron and smuon mean that up to four r-jets are possible in jltlh production 
and efe± production. The branching ratios for three and four r-jets in e\e{ and fiip-i 
production are greatly diminished after cuts. 

The question now arises as to the observability of these modes at the Tevatron's Run 
II. The cross sections for inclusive r-jet production before cuts are given in Fig. |6]. All the 
SUSY production modes considered in this analysis are included. Events with more than 
four r-jets are not included in the figure due to their extremely low branching ratios. By 
far the dominant decay mode is the 2 r-jets mode. The 3 r-jets and 1 r-jet modes are also 
quite large. 

Of course, the real issue is what the production cross sections are after the cuts have been 
imposed. These are given in Fig. [7|. The graph shows that after cuts, the 1 r-jet mode is 
dominant. The 2 r-jets mode is of the same order of magnitude and the 3 r-jet mode is not 
unappreciable. For A = 35TeV, the three r-jets rate is 4.7 fb. For an integrated luminosity 
of 2fb _1 (approximately what is expected initially during Run II), this corresponds to ~ 9 
observable events. For 30 fb -1 , the number of observable events is 141. The 2 r-jets cross 
sections of 48.1 fb gives ~ 96 events for 2fb~ 1 of data and ~ 1440 for 30 fb -1 of data. 
As A increases, the numbers are smaller due to the smaller SUSY production rate. For 
A = 50 TeV, the expected number of events for three r-jets is about 2 for 2 fb -1 of data and 
30 for 30 fb -1 of data. The expected number of 2 r-jets events is 16 and 248 for 2fb _1 and 
30 fb -1 of data, respectively. For A = 70 TeV, the expected number of events for two r-jets 
is 2 and 28 for 2fb _1 and 30 fb -1 , respectively. 

The branching ratios for some of the more interesting individual modes in combined 
SUSY production are given in Table [V]. The electrons and muons are typically the soft to 
pass the cuts. Thus requiring an e or ji to enhance the signal over background probably will 
be of little help. 



V. STAU NLSP LINE WITH N = 3 

We now consider a case where the ordering of the sparticles masses is quite different from 
the previous case. The parameters taken here are n = 3, tan/5 = 15 and M/A = 20. We 
vary A from 25 to 55 TeV. The masses for these parameters are given in Fig. |8|. We see that 
the ordering of the masses is of type 2: M x o > M x ± > ma > M x o > m^,^ > M fl . This case 
is more complicated than the previous one due to the shifting of the sneutrino masses below 
that of xt an d X*2- As a consequence, there are many decay modes for xt an d X2 over the 
parameter space considered here. Moreover, the lightest selectron and smuon masses are 
always below that of the lightest neutralino. The result of all this is that the decay chains 
will generally be quite involved with many steps for the values of A considered here. 



The branching ratios for the sparticles are given in Table [VI|. Since the masses of the 
lightest selectron and the lightest smuon are always below that of the lightest neutralino, 
there are three decay modes available for the values of A considered: Xi ~^ ^iT? Xi fiif^ 
and Xi ~ eie. The decay to the stau is the dominant decay mode especially at low values 
of A. 

There are many potential decay modes for the chargino with these values of the parame- 
ters. Since the sneutrinos are now less massive than the chargino, these provide three decay 
modes that were not present in the previous case: xt ~ > v t t, xt ~ > ^VA* and xt ~^ & e e. 
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For the entire range of parameters considered, these decays to the sneutrinos are always 
present as well as the decay Xi ~^ ^\ v t which is the dominant decay mode for A less than 
about 50TeV. As A increases, the mass difference between the lightest chargino and the 
lightest neutralino increases and the decay xf ~ * Xi W becomes kinematically allowed. At 
A = 55TeV, it is as important a decay mode as Xi T\V T . One other distinguishing 
characteristic of this case from the last one is that the heavier selectron and smuon have 
masses below that of xt an d Xt Thus the decays xt ~* t 2 v t , xf — ► and Xi ~ > Zi v e 
are available and their branching ratios are small, but not unimportant at large A. 

For the second lightest neutralino, there are up to eleven main decay modes. For low 
values of A the dominant decay mode is X2 ~ * ^i r - The other slepton decay modes x\ ~^ 
e\e and x\ ~> fiif^ are a l so important. As A increases, the decay X2 ~> Xih becomes 
kinematically allowed and rapidly becomes the dominant decay. The decay modes to the 
sneutrinos also become more important. 

The Et distribution of the leading r-jet for A = 25 TeV is given in Fig. |1(] and the Et 
distribution of the secondary r-jet is given in Fig. [11]. The distribution for the leading r-jet 
is quite similar to the previous case, but the secondary r-jet spectrum is softer due to the 
decrease in the direct production of r leptons from chargino and neutralino decays and more 
of the r leptons coming from further down the decay chain. The Ifix distribution is given in 



Fig. 12 



We now consider the details of the various final state possiblities. Table |VII] gives the 
inclusive branching ratios for different numbers of r-jets for A = 25 TeV. In principle, up to 
six r leptons can be produced in XiXi production, but the five and six r lepton branching 
ratios are small. The most important mode before cuts is the two r-jet mode at 40.4%, 
but the one and three r-jets modes are also appreciable. After implementing the cuts, the 
branching ratios are greatly decreased: the two r-jets branching ratio becomes only 11.4% 
and the one r-jet branching ratio becomes 23.4%. The three r-jet branching ratio becomes 
essentially negligible. The situation changes as A increases. Table |V111| gives the inclusive 
r-jet branching ratios for A = 40 TeV. The branching ratio for greater numbers of r-jets are 
now larger. This is due to the decrease in the branching ratio for xf ~~ * T\v T from which 
one can get only one r-jet from the chargino and the increase in xf ~^ XiW which can 
give three r-jets and xf ~ * which can also give three r-jets. The two r-jets mode in 
X1X1 production is still dominant at 34.8%, but now the three r-jets mode is appreciable 
at 28.6%. After cuts, the three r-jets rate drops to 5% and the one r-jet mode becomes 
dominant at 30.1%. 

Turning now to X2 xf production, we see that at low A, there is the potential to produce 
up to five r-jets (three from xf ~ * v t t with the subsequent decays v T — > Xi v t and Xi ~^ T\T 
and two r-jets from X2 v T v T ), but the branching ratios for more than three r-jets are 
rather small. As usual the dominant decay mode is to two r-jets with a before cuts branching 



ratio of 40.8%, but the three r-jets branching ratio is large at 29.1% as shown in Table |VT1 



After cuts these fall to 16.8% and 3.3%, respectively. The one r-jet mode becomes dominant 
with a branching ratio of 24.5%. Table |V111| gives the results for A = 40 TeV. The four r-jets 
mode has a substantial decay rate before cuts, but this mode becomes negligible after cuts 
due to the softness of the fourth r-jet due to its production further down the decay chain. 
On the other hand, the three r-jets branching ratio is now higher at 6.7%. 

Slepton production for this case is largely the same as in the previous case. With f\ — > r G 
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being essentially the only decay mode for the lighter stau, the r-jet branching ratios before 
cuts are completely dictated by the hadronic branching ratio for the r lepton. For fifjii 
and efei production, up to four r-jets can be produced, but after cuts the rates for three 
and four r-jets are greatly reduced. 

Putting all the pieces together, we can now answer the question as to the probability 
of observing these events at Tevatron's Run II and Run III. Fig. [13| shows the branching 
ratios for the inclusive r-jet modes before cuts for all the considered SUSY production modes 
combined. The two r-jets mode is the mode with the largest a ■ BR, but the production 
rates for one and three r-jets are close to this. After including cuts, the one r-jet mode is 
dominant and the two r-jets mode is respectably high as seen in Fig. [TJ]. For A = 25 TeV, 
the three r-jet rate is 5.4 fb giving ~ 11 for 2fb _1 of data and ~ 162 events for 30 fb -1 of 
data. The two r-jet rate of 39.3 fb gives ~ 79 and ~ 1179 events, respectively. For the 
higher A value of 50 TeV, the rates are cut down significantly. The two r-jets rate of 1.0 fb 
gives ~ 2 and ~ 30 events for 2fb _1 and 30 fb -1 of data, respectively. 



VI. CO-NLSP CASE 

The co-NLSP case |J refers to when the mass difference between e\ (and and f\ 
is less than the mass of the r lepton. When this is the case, the three-body decay mode 
e.\ — > erfi is not kinematically allowed and the main decay mode for the selectron is the 
two-body mode e\ —>■ e G. The parameters for the example of this case that is considered 
here are n = 3, tan/3 = 3 and M/A = 3. A is varied from 25 to 65 TeV. The masses of the 



sparticles of interest are given in Fig. 15. We see that the oredering of the masses here is 
a special case of type 2: M x o > M x ± > rriy > M x o > ~ m fl . With this ordering of 

the masses there are typically many decays mode of the sparticles to consider. The cross 
sections for the SUSY production modes are given in Fig. [16|. Due to the rapid increase 
in the sparticle masses (especially the gaugino masses) as A is increased, the cross sections 
tend to decrease fairly rapidly. 



The branching ratios for the sparticles of interest are given in Table |XJ. Since the lighter 
selectron and the lighter smuon have about the same mass as the stau, the branching ratios 
for Xi —> Tit, Xi ~^ fiif 1 an d Xi ~^ &i e are nearly equal. The decay to the stau is slightly 
favored. 

The chargino's decays strongly depend on the value of A. For values of A that are 30 TeV 
and below, the dominant decay mode of the chargino is xf ~^ T\v T . As A is increased above 
this, the chief decay modes are the decays to the sneutrinos and the decay xf ~* XiW which 
tends to dominate when kinematically allowed. As A increases above 30 TeV, the masses of 
the heavier sleptons (f 2 , J12 and e 2 ) fall below the mass of the chargino and so the decays to 
these heavier sleptons are allowed as well. 

For the second lightest neutralino, there are again up to 11 main decay modes. At low 
values of A, the decays to the lighter sleptons are dominant with x% having a slight 

edge over the other two slepton decays. As A increases the decays to the sneutrinos gradually 
become more important. In addition the decay X2 ~* X?h an d the decays to the heavier 
sleptons become kinematically allowed and dominate over the other decays. 

In chargino pair production at low values of A, two r leptons are always produced 
because essentially the only decay mode for the chargino is xf ~^ ^\ v t while the stau decays 
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via fi — > r G. On the other hand, the classic three r signature for X2 Xi production will 
be diminished since X2 ~~ > ?iT, X2 ~~ > Ai/- 4 an d X2 ~" > are a ^ roughly equal. Since the 
subsequent decays of the selectron and smuon to the gravitino produces no r leptons (unlike 
the three-body decay modes t\ — > erfi and [i\ — > \irf\ that were dominant in the other 
two cases), there will tend to be a depletion in r-jets here relative to the previous type 2 
case which didn't satisfy the co-NLSP condition. For larger values of A, the situation is 
more complicated, but the decay will frequently involve the lightest neutralino. The lightest 
neutralino in turn tends to decay to f\, Ji\ and t\ roughly equally. Thus there is again a 
relative depletion in events with r-jets. 

The Et distribution for the leading r-jet when A = 25 TeV is given in Fig. [17|. The 
distribution of the secondary r-jet is given in Fig. [Tj|. Qualitatively, these are much the 
same as in the previous cases. At A = 25 TeV, the decay chains are relatively short and the 
r-jets tend to be quite hard. The JpT distribution is given in Fig. |19[ 

We now consider the details of the various final state possibilities. Table [TX| gives the 
inclusive branching ratios for different numbers of r-jets for A = 25 TeV. With xf ~^ 
f\V T being the only decay mode here, X1X1 production produces two r leptons. Thus the 
probability for r-jets before cuts is dictated by the hadronic branching ratio of the r lepton. 
Including cuts diminishes the number of events with a given number of r-jets. For example, 
the branching ratio for 2 r-jets falls from 42% to 10.5%. When A is increased, the situation 
changes dramatically. Table [X] gives the inclusive branching ratios for A = 40 TeV. The 
possibility exists to create many r-jets, but the probability for creating more than three is 
low. In addition, the probability for producing no r-jets is high at ~ 35%. After cuts, the 
only appreciable modes are the two r-jets mode at 10% and the one r-jet mode at 18%. 

Turning now to X2 Xi production, we see that at low A, there is the potential to produce 
up to three r-jets. The rates are diminished by the strong presence of X2 ~* an d 
X2 ~ * A^i/^ however, and the rate for no r-jets is high at ~ 25%. After cuts, the two r-jets 
branching ratio is only about 7% and the one r-jet rate is about 25%. For A = 40 TeV, 
the potential exists to create many more r-jets, but the one and two r-jets modes remain 
dominant with after cuts branching ratios of 19% and 11%, respectively. 

For f^rf , the situation is pretty much the same as it is in the previous cases considered. 
With f\ — f r G being the only decay mode of the lighter stau, the probability for a given 
number of r-jets is completely dictated by the hadronic branching ratio of the r lepton. The 
branching ratios after cuts are largely dictated by the mass of the stau. 

We now consider the possibility of observing these events at the Tevatron's Run II and 
TeV33. Fig. |20| shows the combined production rates for the inclusive r-jet modes before cuts 
for all the SUSY production modes considered. We do not include the cross sections for more 
than three r-jets as these are prohibitively small. In sharp contrast to the previous cases, 
the most typical situation is that no r-jets are produced. For low values of A (A < 40 TeV), 
however, the production rates for one and two r-jets are comparable. The results after cuts 
are shown in Fig. [2l|. The one r-jet mode is dominant and the two r-jet mode is respectably 
high. For A = 25 TeV, the two r-jets rate is about 28 fb which gives about 56 events for 
2 fb -1 of data and 840 events for 30 fb -1 of data. For A = 40 TeV, we have a lower production 
rate of 1.8 fb. For 2fb _1 of data this corresponds to ~ 4 events, while 30 fb -1 of data gives 
about 54 events. 

The branching ratios for some of the more important individual modes are given in 
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Table [XIl . Unlike the previous cases considered, there is the potentialthat modes with 
specific numbers of charged leptons could be important. We see from the table that for 
A = 25TeV, the rate for an electron and a r-jet is 12.7fb after cuts. For electrons and 
muons combined, this is 25.4 fb. For an integrated luminosity of 2fb _1 , this corresponds to 
about 50 events, while for an integrated luminosity of 30 fb -1 , this corresponds to about 762 
events. A better signal is the 2 e + r-jet and 2 fi + r-jet signals. The combined cross section 
for this is 13 fb. For 2fb~ 1 of data, this corresponds to 26 events. For 30 fb -1 of data, this 
corresponds to about 390 events. 



VII. CONCLUSION 



We have considered the phenomenology of GMSB models where the lighter stau is the 
NLSP and decays promptly. We have looked at a wide range of the GMSB parameter space. 
Typically, in those regions of the parameter space where SUSY production can occur at 
an observable rate, the dominant SUSY production modes are xf Xi an d X2X11 although 
slepton pair production can be significant in those regions of the parameter space where xt 
and X2 are too massive to be readily produced. The decays of the SUSY particles lead to 
events containing two or three r leptons plus large missing transverse energy. Searching for 
the r lepton signals by the hadronic decays of the r leptons to thin jets is complicated by 
the fact that, while primary r-jets can have high Et, the secondary r-jets tend to be rather 
soft. As a result, many of the r-jets tend to be eliminated by the cuts. We've shown that the 
most promising channel is the two r-jets mode, while the three r-jets mode can be important 
at the higher integrated luminosity at Run III. This r-jet is degraded in the co-NLSP case 
where the lighter charged sleptons are nearly degenerate in mass. The missing transverse 
energy associated with the events is large providing a good trigger for these events. Good r 
identification will be extremely important to detect the signal. 
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TABLES 



TART T? T 


Branching ratios of the 


sparticles of interest for the parameters n = 2, 


tan (3 = 15 


ana _m /a — o. 




of the Ai are obtained by replacing the 


o Willi a (i 


in the e~i decays. 










A (TeV) 








F)prav TVTnrlp 


35 


40 


50 


60 


70 


80 


85 


— ► Tl f-y 

Al ' 1 


1 


1 


0.6787 


0.5192 


0.4440 


0.3996 


0.3833 


Y ± _> Y W 
Al Al 




- 


0.3213 


0.4808 


0.5560 


0.6004 


0.6167 


A/9 ^ Ti T 

A2 ' 1 ' 


0.5677 


0.5965 


0.6235 


0.3075 


0.2137 


0.1719 


0.1578 


Yn — > //I /y, 
A2 A* if* 


0.2162 


0.2017 


0.1660 


0.0659 


0.0378 


0.0256 


0.0217 


y9 — > ei e 

A 2 -L ^ 


0.2162 


0.2017 


0.1660 


0.0659 


0.0378 


0.0256 


0.0217 


A2 Al ^ 






0.0446 


0.0318 


0.0251 


0.0219 


0.0207 


v9 — > v9 Vi 

A2 Al 11 








0.5289 


0.6856 


0.7550 


0.7780 


Xi ^ ^i T 


1 


1 


0.8577 


0.6542 


0.5659 


0.5215 


0.5072 


Xi AiM 






0.0711 


0.1729 


0.2170 


0.2392 


0.2464 


Xi he 






0.0711 


0.1729 


0.2170 


0.2392 


0.2464 


h -> Xi e 


1 


1 












ej" — ► e~r~r + 






0.5205 


0.5287 


0.5315 


0.5310 


0.5298 


ej" — > e _ r + r _ 






0.4795 


0.4697 


0.4634 


0.4580 


0.4554 


ei — > e G 








0.0016 


0.0050 


0.0110 


0.0148 



TABLE II. Inclusive tau-jet branching ratios for the dominant production mechansims for the 



parameters n = 2, tan/3 


= 15, A = 


35 TeV and M = 


105 TeV. 






Production Mode 


1 r-jet 


2 r-jets 


3 r-jets 


4 r-jets 


5 r-jets 


xfxi- no cuts 
with cuts 


0.4562 
0.2577 


0.4200 
0.1084 








Xl X2 : no cuts 
with cuts 


0.2408 
0.2558 


0.4434 
0.1567 


0.2723 
0.0259 






T+fr : no cuts 
with cuts 


0.4560 
0.2523 


0.4203 
0.0939 








e\e{: no cuts 
with cuts 


0.1128 
0.2383 


0.3118 
0.0778 


0.3834 
0.0003 


0.1766 
negl. 
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TABLE III. Inclusive tau-jet branching ratios for the dominant production mechansims for 
the parameters n = 2, tan/3 = 15, A = 50TeV and M = 150 TeV. 



Production Mode 


1 r-jet 


2 r-jets 


3 r-jets 


4 r-jets 


5 r-jets 


XiXi- no cuts 
with cuts 


0.3194 
0.3355 


0.4099 
0.1610 


0.1627 
0.0100 


0.0291 
0.0004 


0.0027 
negl. 


Xl X2 : 110 cuts 

with cuts 


0.1969 
0.3234 


0.3961 
0.2238 


0.3060 
0.0472 


0.0626 
0.0014 


0.0044 
negl. 


r+rf : no cuts 
with cuts 


0.4561 
0.3345 


0.4201 
0.1370 








efe~i'- no cuts 
with cuts 


0.1130 
0.3199 


0.3119 
0.1207 


0.3833 
0.0023 


0.1765 
negl. 


- 


TABLE IV. Inclusive tau-jet branching ratios for the dominant production 
the parameters n = 2, tan = 15 A = 70 TeV and M = 210 TeV. 


mechansims for 


Production Mode 


1 r-jet 


2 r-jets 


3 r-jets 


4 r-jets 


5 r-jets 


xfxi- no cuts 
with cuts 


0.2333 
0.3647 


0.3821 
0.2119 


0.2549 
0.0322 


0.0728 
0.0023 


0.0078 
negl. 


Xi Xt no cuts 
with cuts 


0.1574 
0.3419 


0.3455 
0.2334 


0.3259 
0.0569 


0.1210 
0.0071 


0.0202 
0.0008 


f^rf : no cuts 
with cuts 


0.4561 
0.3988 


0.4201 
0.1711 








e\e{: no cuts 
with cuts 


0.1165 
0.3901 


0.3133 
0.1495 


0.3795 
0.0013 


0.1743 
negl. 




TABLE V. Production rates in fb for some of the more interesting 
with and without cuts for the parameters n = 2, tan/3 =15 and M/A = 


final state configurations 
3. 






A = 35 TeV 
no cuts 


cuts 


A: 

no cuts 


= 50 TeV 

cuts 


1 r-jet 

e/fi + 1 r-jet 

1 jet + 1 r-jet 

2 jets + 1 r-jet 
e//i+ 2 jets + 1 r-jet 
2 r-jets 




85.44 
157.3 


51.17 
11.91 
21.58 

40.64 


4.98 
9.18 


3.71 
1.59 
2.36 
2.35 
0.69 
3.91 
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TABLE VI. Branching ratios of some of the sparticles of interest for the parameter set with 
n = 3, tan/3 = 15 and M/A = 20. 



A (TeV) 



Decay Mode 


25 


on 

30 


A C\ 

40 


50 


55 


± 

Xl T X V r 


0.8062 


0.6330 


0.3514 


0.2211 


0.1839 


± 






0.0118 


0.0493 


0.0643 


± 

Xi e 2 ^ e 






0.023 


0.061 


0.074 


± 

Xl "rT 


0.0729 


0.1097 


0.1379 


0.1450 


0.1467 


± 

Xi ^ e e 


0.0604 


0.0954 


0.1253 


0.1342 


0.1365 


Xi -»■ Xi W 




U.U666 




U.194o 


U.18o2 


n 

X2 r lT 


0.5760 


0.5653 


0.2894 


0.1678 


0.1385 


X° ^27" 






U.U14Z 


U.U4oz 


U.Uo4z 


,0 = „ 
X2 e l e 


0.1667 


0.1335 


U.U478 


0.0207 


0.0151 


n — 
X2 e 2e 






0.0240 


0.0507 


0.0602 


X° v T v T 


0.0317 


0.0583 


0.0760 


0.0811 


0.0845 


X% ~+ 


0.0294 


0.0547 


0.0722 


0.07764 


0.0810 








0.0196 


0.0136 


0.0119 


x^-x?h 






0.3127 


0.3961 


0.3983 


Xi 


0.7955 


0.6011 


0.4707 


0.4268 


0.4150 


Xi e x e 


0.1023 


0.1994 


0.2646 


0.2866 


0.2925 


e~ — > e~r~f + 


0.5563 


0.5746 


0.5898 


0.5960 


0.5977 


e~ — > e~r + r~ 


0.4437 


0.4254 


0.4102 


0.4040 


0.4023 



TABLE VII. Inclusive tau-jet branching ratios for the dominant production mechansims for 
the parameters n = 3, tan = 15, A = 25 TeV and M = 500 TeV. 



Production Mode 


1 r-jet 


2 r-jets 


3 r-jets 


4 r-jets 


5 r-jets 


XiXi- no cuts 


0.3597 


0.4040 


0.1110 


0.0306 


0.0033 


with cuts 


0.2336 


0.1138 


0.0077 


0.0003 


negl. 


XiX2 : n0 cuts 


0.2110 


0.4077 


0.2912 


0.0443 


0.0085 


with cuts 


0.2452 


0.1682 


0.0332 


0.0010 


negl. 


f+rf : no cuts 


0.4555 


0.4210 








with cuts 


0.2258 


0.0797 








efe^: no cuts 


0.1128 


0.3121 


0.3833 


0.1764 




with cuts 


0.1998 


0.0779 


0.0063 


0.0002 
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TABLE VIII. Inclusive tau-jet branching ratios for the dominant production mechansims for 
the parameters n = 3, tan = 15 A = 40 TeV and M = 800 TeV. 



Production Mode 


1 r-jet 


2 r-jet 


3 r-jets 


4 r-jets 


5 r-jets 


vtv, ■ no cuts 


0.1890 


0.3484 


0.2863 


0.1155 


0.0223 


with cuts 


0.3011 


0.2123 


0.0500 


0.0057 


0.0003 


Xl X2 : 110 cuts 


0.1457 


0.3342 


0.3301 


0.1395 


0.0236 


with cuts 


0.2991 


0.2343 


0.0672 


0.0082 


0.0008 


r+rf : no cuts 


0.4556 


0.4208 








with cuts 


0.3396 


0.1395 








efe^: no cuts 


0.1128 


0.3121 


0.3833 


0.1765 




with cuts 


0.3213 


0.1270 


0.0048 


0.0001 





TABLE IX. Inclusive r-jet branching ratios for the various production mechanisms for the 



parameters n = 3, tan (3 


= 3, A = 


25 TeV and M = 


75 TeV. 






Production Mode 


1 r-jet 


2 r-jets 


3 r-jets 


4 r-jets 


5 r-jets 


xfxi- no cuts 
with cuts 


0.4562 
0.2514 


0.4200 
0.1053 








Xi Xt no cuts 
with cuts 


0.5088 
0.2480 


0.1514 
0.0666 


0.0935 
0.0153 






T\T^\ no cuts 
with cuts 


0.4560 
0.2427 


0.4203 
0.0891 









TABLE X. Inclusive r-jet branching ratios for the various production mechanisms for the 
parameters n = 3, tan (3 = 3, A = 40 TeV and M = 120 TeV. 



Production Mode 


1 T-jet 


2 r-jets 


3 r-jets 


4 r-jets 


5 r-jets 


XiXi- no cuts 


0.2903 


0.2272 


0.0937 


0.0299 


0.0054 


with cuts 


0.1826 


0.0996 


0.0227 


0.0034 


0.0002 


Xi xl- no cuts 


0.2613 


0.2299 


0.0826 


0.0298 


0.0057 


with cuts 


0.1903 


0.1054 


0.0232 


0.0038 


0.0004 


f+rf : no cuts 


0.4560 


0.4202 








with cuts 


0.3419 


0.1410 
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TABLE XI. Branching ratios of some of the sparticles of interest for the parameter set with 
n = 3, tan/3 = 3 and M /A = 3. 



A (TeV) 



Decay Mode 


25 


OA 

30 


4 n 

40 


50 


60 


65 


± 


1 


0.3416 


0.0225 


C\ C\f\f\C\ 

0.0099 


0.0063 


0.0053 


± 

Xi ^r 




0.2147 


0.142b 


0.1390 


0.1402 


0.1410 


± 

Xi ^ e e 




0.2219 


0.1418 


0.1383 


0.1396 


0.1405 


± T/T7 

Xi -►Xi^ 






0.4514 


0.343b 


0.2781 


0.2552 


± 

Xl 






t~\ nor*/ 1 

0.0326 


0.0765 


0.0984 


0.1055 


xf -»■ e2e 






0.0336 


0.0*72 


n noon 


(J.lUbU 


n 

X2 ~> T lT 


0.3340 


0.3135 


0.1384 


0.0577 


0.0302 


0.0232 


x° g i e 


n QOOQ 


A 0(177 


n i £ 
U.lZOO 


U.U5U0 


U.U252 


A A1 OA 

U.U19U 


X2 ~ »■ ^r^r 


n AA7i 
U.UU ( 1 


n aqa/I 
U.UoU4 


U.USZU 


n nnev 


A 1 AA9 

U.IUDZ 


A 1 AQK 




X2 


n nn^7i 

0.0071 


n nono 

0.0303 


0.0819 


0.0986 


U. 11)51 


/ 1 -i no r 

U.1U85 


X2 - »■ Xl Z 






a m 9/i 
U.U1Z4 


a nncQ 
U.UUoo 


U.UUDU 


U.UUDo 


X2 Xi h 






0.2146 


0.2763 


0.2800 


0.2779 


X2 ^2T 






U.U44 ( 


U.UoDD 


A 1 A/IQ 
U. lU4y 


All AA 
U.11UU 


X° g 2e 






U.U40D 


U.Uo I L 


A 1 A£1 


A 1 1 A1 
U.11U1 


X? n T 


U.ooyo 


U.o4o ( 


U.oo / o 


A QQK9 


A QQKK 
U.OOOO 


U.OOOO 


X? - ► ei e 


U.oZUO 


U.oZol 


n 3^1 i 


n q 


A ^^99 


A '?' : !9' : ! 
U.ooZo 


— X^ 7 " 


U.UU4 ( 












-> Xl 


r\ nn c o 

0.9953 


n nn /I n 

0.9940 


n noon 

0.9889 


n no7c 

0.9875 


0.9870 


n nocn 

0.9869 


Z> T — > fiW 




0.0060 


0.0111 


0.0125 


0.0130 


0.0131 


-> Xl^e 


0.9928 


1 


1 


1 


1 


1 


Ve -> xfe 


0.0072 












^2 X?T 


0.4640 


0.8960 


1 


1 


1 


1 


T2 -> Xf^r 


0.4155 


0.0983 










-T2 -> X° T 


0.1206 


0.0057 










e 2 - X?e 


1 


1 


1 


1 


1 


1 


ei — ► Ge 


1 


1 


1 


1 


1 


1 



TABLE XII. Production rates in fb for some of the more interesting final state configurations 
with and without cuts for the parameters n = 3, tan (3 = 3 and M/A = 3. 





A = 25 TeV 




A = 30 TeV 






no cuts 


cuts 


no cuts 


cuts 


r-jet 




28.84 




5.39 


2 r-jets 


70.57 


23.11 


8.01 


4.49 


e/n & r-jet 


38.20 


12.71 


4.35 


3.82 


26/2^ & r-jet 




6.48 


3.99 


2.75 
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FIGURES 

500 -, 




30 40 50 60 70 80 90 

FIG. 1. Masses for the sparticles of interest for the line defined by n = 2, tan/3 = 15 and 
M/A = 3. M o rj M ± and M Al rs M gl . 
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1000 




30 40 50 60 70 80 90 
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FIG. 2. Cross section for the important SUSY production processes at the Tevatron for the 
line defined by n = 2, tan/3 = 15 and M/A = 3. The X2xt cross section includes production of 
both signs of the chargino. 
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FIG. 3. The Et distributions of the leading r jet for the parameters n = 2, tan/3 = 15, M/A 
= 3 and A = 35TeV. In (a), no cuts have been imposed. In (b), the \r]\ < 1 cut has been imposed. 
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FIG. 4. The Et distributions of the secondary r jet for the parameters n = 2, tan (3 = 15, 
M/A = 3 and A = 35TeV. In (a), no cuts have been imposed. In (b), the \r]\ < 1 cut on r-jets 
has been imposed. 
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FIG. 5. TfiT distribution of the secondary r jet for the parameters n = 2, tan/3 = 15, M/A = 
3 and A = 35TeV. In (a), no cuts have been imposed. In (b), the Et/pt and pseudorapidity cuts 
on the jets and charged leptons have been imposed. 
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FIG. 6. a • BR before cuts for the inclusive r jets modes for the parameters n = 
and M/A = 3. 
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FIG. 7. a ■ BR after cuts for the inclusive r jets modes for the parameters n = 2, tan/3 = 15 
and M/A = 3. 
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FIG. 8. Masses for the sparticles of interest for the line defined by n = 3, tan/3 = 15 and 
M/A = 20. M x o « M x ± and M Al w M gl . 
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FIG. 9. Cross section for the important SUSY production processes at the Tevatron for the 
line defined by n = 3, tan/3 = 15 and M/A = 20. The X2X1 cross section includes production of 
both signs of the chargino. 
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FIG. 10. The Et distributions of the leading r jet for the parameters n = 3, tan (5 = 15, M/A 
20 and A = 25 TeV. In (a), no cuts have been imposed. In (b), the \r]\ < 1 cut has been imposed 
the r-jets. 
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FIG. 11. The Et distributions of the secondary r jet for the parameters n = 3, tan/3 = 15, 
M/A = 20 and A = 20TeV. In (a), no cuts have been imposed. In (b), \rf\ < 1 cut on r-jets has 
been imposed. 
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FIG. 12. IpT distribution of the secondary r jet for the parameters n = 3, tan/3 = 15, M/A 
= 20 and A = 25TeV. In (a), no cuts have been imposed. In (b), the Et/pt and pseudorapidity 
cuts on the jets and charged leptons have been imposed. 
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FIG. 13. a • BR before cuts for the inclusive r jets modes for the parameters n = 3, tan (3 = 15 
and M/A = 20. 
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FIG. 16. The SUSY production cross sections for the co-NLSP example where n = 3, tan = 3 
and M/A = 3. 
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FIG. 17. The Et distributions of the leading r jet for the parameters n = 3, tan f3 = 3, M/A 
= 3 and A = 25TeV. In (a), no cuts have been imposed. In (b), the \rj\ < 1 cut on r-jets has been 
imposed. 
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FIG. 18. The Et distributions of the secondary r jet for the parameters n = 3, tan f3 = 3, 
M/A = 3 and A = 25TeV. In (a), no cuts have been imposed. In (b), the \rj\ < 1 cut on r-jets 
has been imposed. 
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FIG. 19. IpT distribution of the secondary r jet for the parameters n = 3, tan/? = 3, M/A = 
3 and A = 25TeV. In (a), no cuts have been imposed. In (b), the Et/pt and pseudorapidity cuts 
on the jets and charged leptons have been imposed. 
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FIG. 20. a ■ BR before cuts for the inclusive r jets modes for the parameters n = 3, tan (3 = 3 
and M/A = 3. 
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FIG. 21. a ■ BR after cuts for the inclusive r jets modes for the parameters n = 3, tan /? = 3 
and M /A = 3. 
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